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Abstract: 

Normally the gravitationally bound heavy 0 + ions in the Earth's ionosphere are in a 
diffusive equilibrium. However, when energized to superthermal energies of a few eV transverse 
to the geomagnetic field, the combined effects of the downward gravitational and the upward 
electric and mirror forces produce interesting flow patterns in the vertical direction like in a 
pulsating fountain. This flow pattern is studied by means of a particle-in-cell code. 


It is well known that the topside ionosphere is commonly populated with a variety of ion 
species such as H + , 0 + , and H* . The light ions (H + and H* ) readily escape the downward 
gravitational force even when they have a small characteristic energy of about a few tenths of an 
eV. On the other hand, for the heavy 0 + ions the escape requires some sort of energization. For 
example, an 0 + ion escapes when its energy W>W es = m Q + g 0 R% ! r, where m Q + is 0 + 
mass, g Q is the Earth's gravity at its surface, R e is the Earth’s radius, and r is the geocentric 
distance in the topside ionosphere. Using R e =6371 km , g a = 9.8 ms ~ 2 and r = 7350 km, 
the escape energy W a = 7.5 eV . For W < W es the 0 + ions can be transported to large altitudes, 
but they eventually fall down. If the energization occurs in the degree of freedom parallel to the 
magnetic field, the falling ions are lost into the ionosphere. On the other hand, when the heating 
occurs in the degree of freedom transverse to the geomagnetic field, the 0 + ions are trapped in 
the flux tube and undergo a bounce motion under the influence of the downward gravitational 
force and the upward mirror force associated with the downward gradient in the geomagnetic 
field 5 . The upward polarization electric field associated with the density gradient also plays a 
critical role in the bounce motion. 



We study vertical plasma flow in geomagnetic open flux tubes, which are typical of the 
polar ionosphere. The topside ionospheric base is assumed to be at r -r b =7350 km . We 
simulate the flux tube from r = r b to a geocentric distance r - = 3 1,350 km . In order to 

study the plasma flow, we employ a particle code [1], We inject macroparticles corresponding to 
0 + and IT 1- ions into the flux tube at the boundary at r — r b and subsequently follow their 
dynamics under the influence of the gravity, polarization electric field and the mirror force 
associated with the gradient in the geomagnetic field. The injected particles are chosen from 
Maxwellian velocity distributions with a zero drift for 0 + and a drift of 2V A for H + , where V A 
is the thermal velocity corresponding to the temperature T a of H* ions at the base. We assume 
that 0 + and electron temperatures at r - r b are also T a . When a steady-state polar wind flow is 
set up in the flux tube, we initiate the transverse heating event for 0 + ions at geocentric distances 
r >8311 km . The heating at such altitudes is expected to be caused by wave-particle interactions 
and it is incorporated in a random fashion as described by Brown et al [2]. We present here the 
flow pattern of 0 + ions. The H* ions do not show any interesting feature except for their 
supersonic escape. Since H* ions have extremely low concentration compared to that of 0 + and 
BT 1 *, they are not included in the model. 

We show the bounce motion of the heated ions and its temporal evolution through a 
sequence of phase space plots after a brief heating event of 5 minutes at an average rate of about 
0.04 kT 0 / s . When kT a = 0.3 eV , this leads to a maximum energization of about 15 eV. 
Therefore, some of the ions, having energies W>W e3 escape while those with W < W es are 
trapped. Figs. 2a, 2b and 2c show the initial diffusive state of the 0 + . Note that the left, middle 
and right panels in Fig. 1 are the phase space plots in r-V n , r-V L and r-W planes, where r 
is the geocentric vertical distance, V\\ and Vj_ are the 0 + velocity components parallel and 
perpendicular to the geomagnetic field line, respectively, and W = (Pi I 2 + P±)> that is, it is 
the total kinetic energy. Each dot in the figures represent an 0 + ion. Figures 2a, 2b and 2c show 
that the bulk of 0 + ions are generally confined below an altitude r = 10 4 km. 



The state of the 0 + flow after 6 hours of the brief heating period is shown in Figs 2d, 2e 
and 2f. Note that some ions have reached near the top of the simulated flux tube and have 
Fjj > 0 , implying that they will escape. While those having F[| < 0 are failing. Those ions are 
transported to high altitudes by the mirror force F m x -j m 0 V±B~^ {d B / d r) . The transverse 
heating increases Vj_ and enables the upward flow. But the gravity pulls them downward if they 
do not have sufficiently large energy to escape. 

The falling ions mirror back and expand upward and they appear in clumps; the clumps are 
seen in the form of evolving bands in the r -V\\ plane as seen from Figs 2d, 2g, and 2j. The 
bands give the appearance of a water fountain. The bands can be seen as clusters of ions in the 
r-V j_ plane. V ± decreases with increasing r due to the conservation of the magnetic moment 
p = y mVl / B . In terms of energy, the ion clusters appear in the form of "bananas," which 
simply manifest that the upgoing ions have slightly more kinetic energy than the downflowing 
ions. In both r-V± and r — W planes, the upflowing (F[j > 0) ions are shown in green while 
the downflowing are shown in blue. 

The cluster formation occurs by the modulation in the electric field E due to the falling 
and rising 0 + ions because E depends on the density gradients. This modulation of E also 
causes the difference in the kinetic energies of the upgoing and downgoing ions in any given 
cluster. 

Acknowledgment: This work was performed under the grant NAGW-2903. 

References 

[1] G. R. Wilson, C. W. Ho, J. L. Horwitz, N. Singh, and T. E. Moore, Geophys. Res. Lett., 17, 
263, 1990. 

[2] D. G. Brown, G. R. Wilson, J. L. Horwitz, and D. Gallagher, Geophys. Res. Lett., 18, 1841, 


1991. 



Figure Captions 


Fig. 1 . Phase space plots in r-V^ (left panels), f~V L (middle panels), and r-W (right 
panels) at several times are shown. Note that the velocity components are normalized with 
respect to the thermal velocity of H" 1 " ions which is about 5.4 km/s for kT 0 = 03 eV and 
energy is normalized in terms of kT 0 . In the middle and right panels, blue represents falling O' 1 " 
ions, while green represents upflowing ones. 



